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ABSTRACT 

Giant arcs in strong lensing galaxy clusters can provide a purely geometric determination of cosmological 
parameters, such as the dark energy density and equation of state. We investigate sources of noise in cosmog- 
raphy with giant arcs, focusing in particular on errors induced by density fluctuations along the line-of-sight, 
and errors caused by modeling uncertainties. We estimate parameter errors in two independent ways, first by 
developing a Fisher matrix formalism for strong lensing parameters, and next by directly ray-tracing through 
N-body simulations using a multi-plane lensing code. We show that for reasonable power spectra, density 
fluctuations from large-scale structure produce > 100% errors in cosmological parameters derived from any 
single sightline, precluding the use of individual clusters or "golden lenses" to derive accurate cosmological 
constraints. Modeling uncertainties similarly can lead to large errors, and we show that the use of parametrized 
mass models in fitting strong lensing clusters can significantly bias the inferred cosmological parameters. We 
lastly speculate on means by which these errors may be corrected. 
Subject headings: gravitational lensing — dark matter — cosmological parameters 



1. INTRODUCTION 

There is strong observational evidence supporting the flat 
ACDM cosmological model, in which the energy budget of 
the uni verse is d ominated by a poorly understood dark sec- 
tor (e.g Bennett et al. 20031: iRiess et al.l[2004l) . An important 
question is whether dark energy, which comprises at present 
£2a w 0.7 of the critical density, is truly a cosmological con- 
stant or instead a dynamical field with evolving energy den- 
sity. Considerable effort has been directed towards develop- 
ing means of measuring the redshift evolution of dark energy, 
and one of the most promising methods involves gravitational 
lensing tomography. The basic idea is to exploit the fact that 
the lensing efficiency for a lens at distance D\ and source at 
D s is proportional to DiDi s /D s . By measuring the variation 
of lensing strength with source redshift, one can measure dis- 
tance ratios as a function of redshift, and thereby constrain the 
background cosmology. Tomographic weak lensing is also 
sensitive to the growth of structure, providing an additional 
handle on dark energy (e.g. Hu 2002), although it may be pos- 
sible to separate the effects of growth and geometry in weak 
lensing tomography (tTain & T avlor 2003; Zhang et alJl2003l) 
if a purely geometric determination of the cosmology is de- 
sired. 

This met hod may also be applied in the strong l ensing 
regime (e.g. lBlandford & Naravanl 19921) . iGolse et al.1 (120021) 
discuss in detail the cosmographic prospects for a sample 
of galaxy clust ers, each with mult iple giant arcs at different 
redshifts, and Souc ailet alJ (J2004) apply this method to the 
galaxy cluster Abell 2218 to derive preliminary cosmologi- 
cal constraints. Clu ster surveys, suc h as the MACS survey 
dEbeling etal.1120011) . EDisCS sur vey (iGonzalez et al]l2002l) . 
or SDSS cluster lensing survey (Hennawi et al. 2004|), and 
wide-area imaging surveys like the CFHT Legacy Survey 2 

1 Hubble Fellow 

2 http://www.cfht.hawaii.edu/Science/CFHLS 



or the RCS-2 survey are expected to find (and are finding!) 
large numbers of arc-bearing galaxy clusters. Frequently, 
deep imaging reveals multiple arcs in the sa me cluster - a 
spectacular example is provided by Abell 1689 ( B enftez et alJ 
2002a) 3 . With this imminent explosion in strong lensing data, 
it is now appropriate to consider realistic sources of noise in 
lensing cosmography, with an aim towards determining opti- 
mal observational strategies. For example, is it better to de- 
vote large amounts of telescope time to a few clusters, in the 
hopes of finding dozens of arcs per cluster which completely 
constrain the mass model, or is it better to survey a large num- 
ber of clusters and study only a few arcs per cluster? 

In this paper, we study two principal sources of noise in 
strong lensing cosmography. First, we consider how uncer- 
tainties in the mass modeling of the lensing cluster translate 
into errors in the derived cosmological parameters. Next, we 
investigate the effects of random density fluctuations in the 
line of sight uncorrected with the main lensing cluster. We 
find that density fluctuations can lead to large errors in the in- 
ferred cosmology along any given sightline. This requires the 
use of multiple sightlines; that is many lensing clusters must 
be studied in order to suppress the noise associated with large- 
scale structure. With limited telescope time, this may also 
lead to problems: the limited constraints provided by small 
numbers of arcs per cluster require the use of parametric mass 
models in fitting the lens data, and we show that such para- 
metric models can produce biases in the inferred cosmologi- 
cal parameters. 

2. SOURCES OF NOISE 

It is perhaps worth emphasizing the well-known point that 
determining the properties of dark energy using gravitational 
lensing tomography is quite difficult, because variation of 
dark energy parameters produces small effects on the lensing 
observables. For tomographic lensing, the basic observable is 

3 see also http://hubblesite.Org/newscenter/newsdesk/archive/releases/2003/01/image/a 
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FIG. 1. — Ratio of lensing efficiency, G(z)/G(\), as a function of red- 
shift. The curves correspond to dark energy equation of state parameter 
w = —0.8, —1, —1.2 from top to bottom. 



the variation of lensing efficiency G D(D( S /D S with source 
redshift. We do not measure the absolute lensing efficiency, 
only the ratio of efficiencies at different redshifts. In figure [2 
we plot the ratio G(z)/G{\) for lens redshift z/ = 0.3, as a 
function of source redshift z for several values of w. As can 
be seen, in order to produce interesting constraints on w we 
must be able to measure the lensing efficiency at the percent 
level, which is quite challenging. Weak lensing surveys aim 
to overcome this challenge using brute force, by surveying 
large fractions of the sky and measuring ellipticities of huge 
numbers of galaxies. Because strong lensing events are rare, 
it would be desirable to focus instead on a small number of 
objects and study them well. The fact that the signal is so 
small, however, means that strong lensing cosmography will 
be subject to a wide variety of error sources. 

2.1. Uncertainties in mass modeling 

Uncertainties in mass modeling have long impeded the cos- 
mographic applications of strong lenses. For example, the 
main systematic uncertainty in determination of the Hubble 
constant Hq using lens time del ays continues to b e the radial 
density profile of the lens (e.g. Kochanek 2002). Similarly, 
uncertainties in the mass modeling of lensing clusters also in- 
duce errors in the inferred cosmology. The reason is easy to 
understand - variations in the dark energy evolution, which 
change E cr ; t (z), can be compensated by adjusting the model 
density profile T.(R) to hold fixed the observables, which de- 
pend only upon the ratio E/E cr ; t . 

In the appendix, we calculate the parameter errors expected 
given a lens model and a specified level of observational er- 
rors. As shown in AppendixlAl we can write the Fisher matrix 
as a sum over Fisher matrices for individual sources, 



Ftot = / ^Zsi-c"(Zsrc) / d t<F(u,Z src ). 



(1) 



Here, n(z src ) describes the number density of sources as a 
function of redshift z src , and the integral over source posi- 



tion u is restricted to the strong lensing region at each red- 
shift. To illustrate the level of errors expected, let us con- 
sider stron g lensing by tw o classes of models: an ellipsoidal 
NFW (Navarro et alJI 19971) m odel, and the power-law ellip- 
soid model of Barkana ( 1998). We additionally allow param- 
eters for external shear, along with the cosmological param- 
eters (here taken to be only Q.m and w), all of which must 
be simultaneously determined from the image data. Assume 
that there are N multiply imaged sources used to constrain the 
model, drawn from a population followi ng a red shift distribu- 
tion n{z) « z 2 exp(— (z/z w ) 2 ) (iHu & Jainll2003l) over a range 
0.8 < Zsrc < 4, with median redshift z m ed = 1-5. Unless stated 
otherwise, we will assume for the underl ying cosmology a fla t 
ACDM model with WMAP parameters (Be nnett etaLI 2003): 
Q. M = 0.27, D. A = 0.73, D. B h 2 = 0.024, h = 0.72, n=\, 
cr 8 = 0.9. 

We first consider cosmological constraints from power-law 
ellipsoids. The fiducial parameter values we adopt for this 
lens model are: Einstein radius b = 10" (for z src = 4), axis 
ratio q = 0.7, isothermal profile, vanishing core radius, no 
external shear, and CIm = 0.27, w = — 1. For concreteness, 
we take observational errors of a = 0.1", appropriate for 
Hubble Space Telescope imaging. Ignoring the degenera- 
cies between the cosmological parameters (Q.m,w) and the 
lens parameters we would obtain errors On = 0.07JV- 1 / 2 and 
G w = O.lAf -1 / 2 with N arcs. However, marginalizing over the 
lens parameters, these errors explode up to an = 3.6N~ 1 / 2 
anda K = 5N- 1 / 2 . 

Next we consider the NFW model. One qualitative differ- 
ence between this model and the singular isothermal ellip- 
soid is that the NFW profile produces a central image at small 
radii. This increases the dynamic range of the lensing con- 
straints, and helps break degeneracies between the cosmol- 
ogy and the mass profile. Using fiducial NFW parameters of 
M v i r = 5 x 10 14 /i _1 M Q , c v ir = 5, axis ratio q — 0.7 and no 
external shear, we find improved constraints: Go = 0.5N~ X I 2 , 
(7 H = 0.1N~ X I 2 . The improvement can be explained in part by 
the extra constraints from the central image, and in part by the 
smaller number of parameters in the NFW model compared to 
the power-law ellipsoid. 

In reality, central images are rarely detected in strong 
lenses, and are strongly affected by the properties of the cen- 
tral cD galaxy in the cluster. Additionally, we expect real 
clusters to exhibit a broader range of density profiles than 
can be accommodated by the NFW model, due to, for ex- 
ample, the effect of baryons. We therefore expect the errors 
in real clusters to be significantly larger than the errors we 
have computed for NFW ellipsoids. Note that these errors are 
only those caused by propagating the observational uncertain- 
ties; below we discuss additional sources of error. Already, 
though, these error levels would require large numbers of arcs 
in a given cluster for useful constraints to be placed on dark 
energy. 

2.2. Systematic errors in the mass model 

One reason for the large errors found above is that we em- 
ployed models with many parameters, producing degenera- 
cies with the cosmological parameters. This requires large 
numbers of arcs in order to constrain the model parameters 
sufficiently. Had we taken a less flexible model, we would 
not have found as large errors on the cosmological param- 
eters. Unfortunately, real clusters are expected to exhibit a 
broad range of behavior, as found in N-body simulations. If 
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FIG. 2. — Fits of NFW lenses to softened power-law ellipsoidal models. 
Each point corresponds to the best-fit values of Q.m and w derived by mod- 
eling sets of 50 randomly generated arcs at 5 different redshifts. The input 
cosmology had &m = 0.27 and w= —1, as depicted by the filled circle. 



the assumed mass model does not fully capture the range of 
behavior of the real density profile, then the derived cosmo- 
logical parameters will be biased; that is, the fit can be im- 
proved by adjusting the cosmology to make the convergence 
k(R,Zs) a better fit to the density profile than variations of the 
mass model T.(R) can accomplish alone. 

This bias is illustrated in figure [2] The figure shows re- 
sults of a Monte Carlo calculation, in which artificial arcs 
are generated using ellipsoidal NFW profiles, and fit ted us - 
ing softened power-law ellipsoidal models (Bar kana|ll998l) . 
which can mimic a wide range of radial profiles over the rel- 
evant radii of interest for strong lensing. A CMB prior on the 
D.M — w plane has also been imposed: the distance to the last 
scattering surface at z = 1089 is requir ed to be d& = 14 ±0.2 
Gpc for h = 0.71 (Bennett et afl l2003l) : the degenerate direc- 
tion is along decreasing w (i.e. more negative) as Q.m is in- 
creased. For each realization, we generate 10 arcs each at 
redshifts of z src = 0.8, 1.6,2.4,3.2 and 4 using the same NFW 
model, and fit all 50 arcs using the power-law ellipsoid model, 
allowing Q.m and w to vary in the fit. We plot only the best- 
fitting Q.m and w for each realization; although the confidence 
region for each realization is large, we can detect parameter 
biases by examining the best-fitting models. Clearly, the cos- 
mological parameters derived from multiple arcs can be bi- 
ased at significant levels (in this example, Aw ~ 1). The sign 
and direction of the bias depend upon the profile shapes for 
the cluster and the model, and could be different from clus- 
ter to cluster if there is significant variation among cluster 
profiles. The bias can be suppressed by using more flexible 
models in the fits. Unfortunately, as discussed above, extra 
degrees of freedom in the mass model come at the expense of 
weaker constraints on cosmology. 

2.3. Projection noise 

In addition to the modeling uncertainties discussed above, 
another source of noise in strong lensing cosmography is 



large-scale structure. Fluctuations in the density along the 
line of sight, if uncorrected, produce errors in the derived 
cosmological parameters. Again, this effect is easy to under- 
stand - density fluctuations at redshifts different than the pri- 
mary lens redshift have different lensing efficiencies than the 
primary lens, producing effects varying with source redshift 
which appear similar to variations in the background cosmol- 
ogy. Naively, one might expect the parameter errors to be 
small - the density fluctuations are at the < 1% level, and ge- 
ometric factors suppress the effects of fluctuations at redshifts 
very different from the cluster redshift. 

In the appendix, we compute the errors expected in the cos- 
mological parameters caused by density fluctuations along the 
line of sigh t. We use the same fiducial lens models as in sec- 
tion £12 2.11 and estimate the line-of-sight density fluctuations 
by pro jecting the matter power spectrum using the Limber 
Jl95l ' approximation, 



si? = 



Z2 H 
, ~H 



dk A 2 (k,z) 
k krn 



Wi(krG), (2) 



with T.Q = cH /4nG, r H = c/H (h A 2 = k i P(k)/2n 2 , and 
W2M = 2Ji(x)/x, with a smoothing scale 9 = 30". We 
smooth the power spectrum in order to select only those 
fluctuations that are coherent over the strong lensing region; 
smaller scale fluctuations could be suppressed by observing 
large numbers of arcs in the cluster, and by modeling of indi- 
vidual galaxies as substructur e in the mass map. We us e the 
nonlinear power spectrum of Peacock & Dodds ( 1996) and 
the approximate transfer functions of Eisenstein & Hu ( 1999), 
for a flat ACDM cosmology with WMAP parameters. We 
used 10 lens planes to describe the line-of-sight density fluc- 
tuations, assuming that the shear variance equals the conver- 
gence variance (and assuming that the shear is randomly ori- 
ented in each plane). When we include these line-of-sight 
fluctuations and single-plane lens modeling, as described in 
the appendix, we obtain errors of an ~ 0.34 and ff„. ~ 0.32 
for the NFW model, and On — 1 and a w ~ 2 for the singular 
power-law ellipsoid model. Note that these errors do not di- 
minish when large numbers of arcs are observed per cluster; 
the projection errors are suppressed only by viewing multiple 
independent sightlines. 

3. COMBINED ERROR FORECAST 

In this section, we combine all of the above sources of noise 
to compute the expected error budget for strong lensing cos- 
mography with giant arcs. We r ay-trace thro ugh N-body sim- 
ulations using code described in DalalZinjSxM), modified 
slightly to perform multi-plane lensing, and construct mock 
catalogues of giant arcs. Then we fit the resulting arc data 
using ellipsoidal NFW mass models, optimizing both the lens 
parameters and the cosmological parameters Q.m and w. One 
advantage of using N-body simulations is that we can account 
for the effects of the highly skewed distribution of density 
fluctuations in the nonlinear regime; in principle this non- 
Gaussianity could lead to additional biases in the cosmolog- 
ical parameters which would not be revealed by the Limber 
approximation used above. 

The N-body si mulati on we use is described in 

f imbsganss etafl J2003I) . Briefly, the TPM code 
ode & Ostriker 2003) was used to simulate a box with side 
length L = 320/j~ 1 Mpc in a flat ACDM cosmology with 
Q M = 0.3, £2 A = 0.7, h = 0.7, n s = 1 and cr 8 = 0.95; this is 
broadly consistent with WMAP parameters iBennett et alJ 
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FIG. 3. — (Left) A simulated cluster at z = 0.29 used to generate arcs. The 
cluster has a virial mass of M = 10 14 5 h~ l M@ and c v j, = 7. The side length is 
2/r 1 Mpc comoving. (Right) The z slc = 5 critical lines for this cluster. The 
black curve corresponds to the critical lines for no intervening lens planes, 
and the colored curves are the critical lines for 10 different realizations of the 
intervening matter. 
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120031) . A total of N ~ 1024 3 particles were used in the 
simulation. We then sliced through the simulation volume to 
produce lens planes spaced every 160/i _l Mpc with angular 
size roughly 20', so that the physical side length of the 
projected planes telescopes with increasing distance from 
the observer at z = 0. A total of 243 pairs of planes were 
produced for each of 19 redshifts between z — and 6.37. We 
subtracted the cosmic mean density from each plane and then 
tiled the lightcone by randomly selecting a pair of planes for 
each redshift, for a total of 38 lens planes per realization. 

In order to ensure that strong lensing occurs, we insert an 
extra lens plane containing a massive cluster, along with all 
the surrounding matter within 5/i _1 Mpc. An example of one 
of the clusters used is shown in fig. [5] We have used only 
clusters with regular, relaxed cores, since clusters with mul- 
tiple massive components require more sophisticated model- 
ing than is feasible using automated routines. We first gener- 
ate arcs at multiple source planes using only the cluster lens, 
neglecting the other lens planes tiling the light cone. This 
amounts to assuming that matter is unclustered at redshifts 
different than the cluster redshift. We fit these arcs using the 
projected NFW model, holding fixed the cosmological param- 
eters, to create a starting model from which subsequent fits 
are started. Next, we tile the lightcone by randomly select- 
ing slices as described above, and recompute the ray-trace. 
We again generate sets of strongly lensed sources, and model 
the mock data sets using single-plane mass models, starting 
from the models which best fit the single-plane lens data, op- 
timizing over both lens parameters and the cosmological pa- 
rameters Q.m and w. In general, the lensing properties of the 
clusters remain quite similar with and without the extra lens 
planes. For example, the right panel of fig.|3]shows the criti- 
cal curves for the cluster shown in the left panel, for sources 
at Zsrc = 5, and it is apparent that the critical curves are only 
slightly affected by the LOS density fluctuations. Similarly, 
the total lensing cross sections and optical depth are negligi- 
bly affected by inclusion of the LOS material. 

Nevertheless, considerable errors in the derived cosmolog- 
ical parameters can arise from projection noise as well as de- 
viations of the cluster profiles from the simple NFW form. 
Figure @] shows results from our ray-tracing simulations. As 
is apparent from the figure, quite large, biased errors can 
arise from the combination of effects discussed above. The 
95% confidence regions for the two fitted parameters are 
0.13 < D. M < 0.45 and -3.3 < w < -0.7 with mean values 
of (£Im) —03 and (w) = —1.3, for input values of Q.m = 0.3 
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FIG. 4. — Cosmological parameters derived from fitting sets of 50 strongly 
lensed sources between redshifts z slc = 0.8 to 5, generated by ray-tracing 
through N-body simulations with full lightcone tiling. Each point corre- 
sponds to a realization of the intervening planes and lensed images. An input 
cosmology of Q.m — 0.3, w = — 1 was used to generate the lensed images. 



and w = — 1 . 



4. CLEANING OUT THE NOISE 



The sources of noise discussed above pose serious chal- 
lenges to the use of giant arcs for cosmography, however they 
are not necessarily insurmountable. Modeling uncertainties 
can be removed by observing large numbers of arcs in each 
cluster, so that a nonparametric, unbiased reconstruction of 
the mass profile and cosmology may be performed. This re- 
quires high resolution and deep imaging, along with redshifts 
for the arcs. 

Projection noise appears more recalcitrant. Directly mea- 
suring the intervening mass fluctuations with weak lensing 
(e.g.|javlor 2001) is difficult, in part because the main lens 
cluster must be carefully subtracted, but also because weak 
lensing can only measure the fluctuations on scales larger 
than a few arcminutes, while the strong lensin g region is 
sensit ive to power down to sub-arcminute scales ( Dal ai et all 
2003). Another possible means of cleaning out the projection 
noise would be to reconstruct the density and shear fluctu- 
ations along the line-of-sight using the observed galaxies in 
the field. This is not unreasonable, since we know that most 
of the pow er on the re levant ^arcminute scales is produced 
by halos (Seliak 2000). However, previous attempts at such 
reconstructions have not been entirely successful, with differ- 
ent groups finding quite differ ent shear and convergence over 
the same field (e.g.lBeiute z et alJ2 002bHLewis & Ibatal2001t 
Morts ell et all200lURiess et al!2001l) . 

The main difficulty with predicting the lensing convergence 
and shear fields from the field galaxies is that galaxies and 
mass are not perfectly correlated. Recall that we must be able 
to suppress the projection noise by at least 1-2 orders of mag- 
nitude in order to derive useful cosmological constraints from 
multiple arcs in a single cluster. Unfortunately, the cross- 
correlation between galaxies and mass, r — P gm / {P gs ,P m m) X I 2 , 
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is expected to be less than unity, typically of order r ~ 0.5 
on the relevant scales. 4 This would appear to preclude a suf- 
ficiently accurate reconstruction. Another way to understand 
this is as follows. Reconstructing the convergence and shear 
arising from individual galaxies in the field involves associat- 
ing halos of mass M with galaxies of luminosity L. Unfortu- 
nately, the relation between M and L is not one-to-one. For ex- 
ample, a galaxy like the Milky Way may be the central galaxy 
in halo of mass 1O 12 M , or it may be a member of a group 
an order of magnitud e larger in mass . In th e halo occupation 
distribution model of iRravtsov et alJ fc00 4). there is roughly 
a factor of 3 — 5 uncertainty in the mass of a halo hosting a 
given galaxy, which directly translates into the uncertainty in 
the convergence and shear associated with the galaxy. Even if 
central galaxies may be distinguished from satellite galaxies, 
there is still an unknown scatter in the relat ion between halo 
mass and central ga laxy luminosity ( Abazaiian et al.ll2004t 
Tasitsiomi et al. 2004). In addition, halo triaxiality and sub- 
structure will affect the convergence and shear at the level of 
interest. 

One crude test of the reconstruction method may be possi- 
ble using extant cosmic shear data. Using the observed fore- 
ground galaxies in the field, one can make a prediction for the 
cosmic shear as a function of position, which can be checked 
against the values actually observed. This unfortunately can- 
not test the fidelity of the reconstruction on the required ~ 10" 
scales, since shot noise in the source galaxies typically lim- 
its the resolution of shear maps to arcminute scales. While 
this cannot fully validate the method, a test along these lines 
would be a worthwhile indicator of whether significant ob- 
servational efforts should be devoted towards the use of giant 
arcs to constrain geometry. 

5. SUMMARY 

We have quantified various sources of noise expected in 
strong lensing cosmography. The use of multiple arcs in 
galaxy clusters holds great promise towards constraining the 
properties of dark energy, however significant errors can arise 
from line-of-sight projections and modeling uncertainties. 

4 Note that observationa l estimates (e.g. Hoekstra et al. 2001; Seli ak et a!] 
120041 ISheldon et al. 2004) typically quote r > 1 on small scales because the 



Removing the projection noise to the necessary level appears 
difficult, since reconstruction of the mass density along the 
line of sight is hampered by the fact that mass and light are 
not perfectly correlated. If projection noise cannot be ade- 
quately cleaned out, then we cannot use "golden lenses", but 
must instead average over many strong lensing clusters in or- 
der to derive useful cosmological constraints. To avoid biases 
in the inferred cosmological parameters, highly flexible or 
nonparametric mass models must be used to describe the clus- 
ter. In the absence of any knowledge of cluster density pro- 
files, non-parametric models require large numbers of arcs to 
be observed in each cluster to sufficiently constrain the mass 
model. However, in principle one could use N-body simu- 
lations to impose prior constraints on cluster profiles. Using 
priors from N-body simulations would not be useful for indi- 
vidual systems, because of the large scatter in profiles among 
individual clusters. However, if simulated clusters do provide 
a reasonable description of the population of real galaxy clus- 
ters, then priors derived from simulations may allow useful 
cosmological cons traints to be derived from large ensembles 
of lensing clusters. Meneghett i et al J (120041) make preliminary 
investigations into this approach. It may also be worthwhile 
to investigate whether additional information, such as X-ray 
temperature profiles of lensing clusters, can be used to im- 
prove the cosmological constraints derived from giant arcs. 
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shot noise is subtracted from the galaxy power spectrum P s 



APPENDIX 
A. COVARIANCE MATRIX 

Consider a model for a lens system with « D b s observables O and n p parameters p describing the lens (and possibly cosmology 
as well). Assume that the penalty function may be written x 2 = SO ■ V • SO, where V -1 = (SO SO) describes the covariance of 
the observational errors. For each lensed source, there are additionally m parameters, with m = 2 if only the image positions are 
being modeled, or m = 3 if the fluxes are modeled as well. The mn src source parameters are nuisance parameters, over which we 
marginalize to obtain the n p x n„ covariance matrix and its inverse, the Fisher matrix F. If D = dO/dp, then propagating errors 
we have F = D T ■ V • D. For a single source, the (unmarginalized) Fisher matrix takes the form 

Hp m 

np \ (T ^\ (ai) 

m{ \ B C J 



where the n p x n p block A corresponds to the lens parameters, the m x m block C corresponds to the source parameters, and 
the n p x m block B is the cross term. Marginalizing over the source parameters, the resulting Fisher matrix becomes F = 
A-B C 1 B T . 
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For multiple sources, the full (unmarginalized) Fisher matrix has the form 



/LA; 



BJ 



Bi 
Ci 






B 2 



c 2 





B 3 



c 3 



V 



(A2) 



Here, A,-, B,, and C, are the components of the single-source Fisher matrices (ean. lAlt for source i. We wish to invert F„ to 
obtain the parameter covariance matrix. We can do so by noting that F„ is diagonalized by the transformation 



U-F„-U T 



where 



U: 



/ 1 
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/F tot 
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B 2 C 2 
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and F t , 



= LFj = E/A J -B i .c; 
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From this, it is easy to see that 
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The covariance of the n p lens parameters is given by the upper left n p x n p block F tot . Therefore, the full (many-source) 
marginalized Fisher matrix is simply the sum of the marginalized single-source Fisher matrices, as we might have expected. 

B. LARGE-SCALE STRUCTURE 

Density fluctuations along the line of sight will perturb the lensed images, inducing changes in the derived best-fitting lens 
parameters. iDalal & Kochanek (2002) compute the perturbations to the parameters to first order in the density and shear pertur- 
bations, which is sufficient for our purposes (recall that the rms c onvergence fluctuations are (k -2 ) 1 / 2 ~ 1%). The only minor 
modification here is that multi-plane lens factors must be included (lKeetonl2003t fSchneider 'et al.1 1 9921) . 

Let us parametrize the line-of-sight fluctuations by projecting the density and shear onto additional lens planes. Each plane 
has its own surface density perturbation 5E (measured relative to mean density) and (2-component) tidal shear 8T, so that for N 
additional planes we have n q = 3N variables, which we denote by n q element q. Just as shifts in the lens parameters Ap produce 
shifts in the observables (dO/dp)Ap, the line of sight fluctuations produce shifts (dO/dq)Aq. Minimizing % 2 , this implies that 
the LOS fluctuations induce parameter shifts 

T 



Ap = F,7 1 .( 4^ 



dO 
dp 



d0 A 

— Aq 
dq 



= M Aq. (Bl) 

We are mainly interested in the induced shifts in the lens parameters and cosmological parameters, not the nuisance parameters 
describing the sources. Let us write the first n p elements of the full para meter array p as pi en s, and the m nuisa nce p arameters for 
each of the sources as u. Using the above expression for F u (ean. lA5> . we can write the first n p rows of ean. lBll as a sum over 
independent terms for each source, 

T / -,^,\ T" 



M 



^-1 



( dOj 

\dp 



lens 



* du 



•V- 



dOj 
lq~ 



(B2) 



Given our expression for M, the parameter covariance caused by LOS fluctuations is (ApiApj) = M; a (Aq a Aqp }Mjp . 
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